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Chapter I 

Introduction 


Pitot tubes are used extensively in fluid flow to measure total pressure for velocity calculations. It is well known, 
however, that the presence of the pitot probe in shear flow — even in incompressible flow — alters the nature of the 
flow. The nature of the alteration can be divided into three effects: 

1. Reynolds Number Effect: This effect is important when R D < 200 (Reynolds number based on pitot tube 
diameter). 

2. Wall Effect: This effect is only seen when the pitot tube is stationed within two diameters of a solid boundary. 

3. Velocity GradientyCenterline Displacement Effect : This effect is always present, and is caused by two things: 
(1) the pitot probe deflecting the incoming streamlines and (2) the average pressure across the pitot probe face 
is not equal to the impact pressure measured at the geometric center of the probe. Figure 1.1 illustrates the 
centerline displacement effect. 

Of these three effects, the centerline displacement effect is the most predominant in high speed boundary layer 
flow. It has been the goal of previous researchers and is the goal of this paper to understand the nature of the center- 
line displacement effect in order to compensate for it. 


NOMENCLATURE 

U.S. Customary Units were employed for the experimental measurements in this study, but the International Sys- 
tem of Units (SI) are used to report the results. 

C constant in law-of-the-wall equation 

C f coefficient of friction 

D external diameter of pitot probe 

d internal diameter of pitot probe 

M Mach number 

p s static pressure 

p t total pressure 
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Pj 2 pitot pressure 


R d Reynolds number based on D, (p e u e D)/p e 
T t total temperature 

T s static temperature 

U + dimensionless velocity (u*/u T ) 

u local velocity 

u x frictionless velocity ( ^/x^/p ) 

u* Van Driest generalized velocity 

y distance from the wall 

y + dimensionless distance from the wall (yu^u) 

Greek 

y ratio of specific heats, c p /c y = 1 .4 for air 
A centerline offset 

8 boundary layer thickness 

k Von Karman constant 

p dynamic viscosity 

v kinematic viscosity 

p density 

x shear stress 

Subscripts 

e property evaluated at the boundary layer edge 
wall property evaluated at the wall 
oo property evaluated in the free stream 

0 wind tunnel plenum conditions 

1 property evaluated upstream of shock wave 

2 property evaluated downstream of shock wave 



1.1 LITERATURE REVIEW 


The pitot probe centerline displacement effect was first studied by Young & Maas (1936). For their experiment 
they made multiple pitot traverses behind a symmetrical wing section in a 0.71 x 0.508 m (2.33 x 1.67 ft) open re- 
turn-flow wind tunnel at Cambridge. They ran two airspeeds: 1 82.9 and 230.1 m/s (60 and 75 ft/s) and used six dif- 
ferent-diameter, square-cut pitot probes with an average inside to outside diameter ratio (d/D) of 0.6. From their 
experiments, Young & Maas discovered that their velocity versus wall height curves were offset from the “true 
curves” (their true curves were determined using extremely small pitot probes). They theorized that this offset could 
be viewed as either a AU or a Ay. They chose to view their offset as a Ay because theory predicted that the incoming 
streamlines would bend around the protruding pitot tube and, thus, be offset by a certain Ay. Young & Maas 
noticed that the effective center of their pitot tubes were displaced from the geometric center towards the region of 
higher velocity. Denoting this centerline offset A, 1 they found that A/D (where D is the outside diameter of the pitot 
tube) was relatively independent of D and was equal to 0.18 when d/D equaled 0.6. Young & Maas then performed a 
second series of tests to determine the displacement effect due to varying d/D (inner diameter/outer diameter). Three 
different pitot tubes with d/D = 0.138, 0.35, and 0.73 were used. From these tests Young & Maas created the relation 
that A/D = 0.131+ 0.082(d/D) for any d. 

The next researcher to investigate the centerline displacement effect was F.V. Davies (1952). Davies’ experiments 
were done in a 0.127 x 0.127 m (5 x 5 in.) supersonic wind tunnel with M^ = 2.43. He ran experiments using both a 
flat plate and a cone, and his flow was laminar for both cases. Davies divided his results into three distinct catego- 
ries: displacement of the boundary layer velocity profile (variation in slope), distortion of the velocity profile near 
the free stream, and distortion of the velocity profile near the boundary surface. Unlike Young & Maas, however, 
Davies found that his profiles were displaced in the direction of lower velocity (an effect he later concluded was due 
to the fact that his probes were flattened probes and Young & Maas’s probes were circular). From his research 
Davies found an error of A/D = -0.15. Davies set a lower limit of d/D = 0.2 on the inner to outer diameter ratio to 
limit the velocity profile distortion near the free stream. 

Livesey (1956) conducted experiments in the turbulent boundary layer of Rolls-Royce’s 7 x 5 ft wind tunnel. His 
purpose was to select a probe design that would have negligible centerline displacement effects. For his experiments 
he used a conical nose probe with d/D = 0.7, a conical nose probe with d/D = 0.5 and a hemispherical nose probe 
with d/D = 0.5 (remember that previous researchers had used square-cut pitot probes). Livesey found that the conical 
nose pitot probe produced a A/D of less than 0.02 while the hemispherical nose produced a A/D of 0. 10 (compared 
to Young & Maas’s 0.15). From these results, Livesey concluded that the conical nosed pitot probe might be the best 
probe to use since it appears to have a negligible displacement error. It has, however, since been concluded that 
conical probes have their limitations in supersonic flow (Bryer & Pankhurst, 1971). 

MacMillan performed his well documented and oft-referenced experiments in 1957. His goal was to determine 
both the centerline offset error and the wall effect error and to combine them into one error or error formula if pos- 
sible. MacMillan ran his experiments using laminar, low-speed flow conditions in a pipe and on a flat plate with 
zero pressure gradient. Five circular, square-cut pitot probes with average d/D = 0.6 were traversed across the pipe 
and through the flat plate boundary layer to measure pressure/velocity versus y. At every y-value the Mach versus 
D was plotted for each pitot tube and the resulting points were curve-fitted then extrapolated back to zero to deter- 
mine the Mach at D = 0. These resulting points were plotted as the true curve and each pitot probe error was deter- 
mined from their offset from this true curve, MacMillan found that when y/D > 2 (the probe centerline was located 
greater than two diameter from the wall, see fig, 1,2), the effective center is displaced towards the region of higher 
velocity by a constant centerline offset of A/D = 0.15 (versus A/D = 0.18 for Young & Maas). When y/D < 2 then an 
additional correction — a u/U correction — must be applied. This new correction is solely a function of y/D and 
MacMillan plotted the resulting curve in what has become a well reproduced curve (fig. 1.3). 

P. Davies (1957) conducted subsonic turbulent boundary layer and wake flow using square-cut pitot probes of 
five various outer diameters. His purpose was to compare the errors of these two types of flow conditions (boundary 
layer versus wake). For his boundary layer testing, Davies theorized that the worst error would occur with the probe 
on the wall and thus took data only with each probe touching the wall. Davies found no appreciable error and thus 
concluded that “the absence of such corrections in the present results may be due to the higher response rate of the 
manometer system used.” 


1 Young & Maas denoted centerline offset by "6,” but here “A” is used to avoid confusion with boundary layer thickness also denoted by “6.” 
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The final and most recent pitot probe centerline displacement research was done by Jerry M. Allen (1972) at 
Langley Research Center. In his Technical Note, Allen outlined an experiment very similar to MacMillan’s, but his 
probe survey was conducted in a supersonic, turbulent boundary layer (as opposed to an incompressible, turbulent 
boundary layer). His purpose was to provide experimental data on the pitot-probe displacement in a supersonic tur- 
bulent boundary layer and to investigate how different ratios of D/8 (where 8 is the boundary layer thickness) affect 
the boundary layer profile. Allen used probes with D/8 ranging from 0.018 to 0.691 and ran his experiments at 
M^ = 2 and a unit = 8xl0 6 /m. Unlike MacMillan and P. Davis, Allen’s A/D may be dependent on D and is 
approximately equal to 0.38 at D = 0 (a calculated value). Allen also found that “large” probes distorted the shape of 
the boundary-layer profile — large being defined as D/8 > 0.145. Figure 1.4 recreates figure 1 1 in Allen’s paper and 
neatly summarizes previous research on centerline displacement error. Table I summarizes the conditions and results 
of previous investigators. 


1.2 SCOPE OF WORK 

The purpose of this research was to continue the investigation of the centerline displacement error for supersonic 
turbulent boundary layers. Data were collected in three NASA Lewis Research Center Supersonic Wind Tunnels 
(SWT’s): the 10 x 10 ft SWT, the 1 x 1 ft SWT and the 15 x 15 cm SWT. Table II summarizes the Mach num- 
bers, Reynolds numbers, D/8, and boundary layers of the present experiment. Included in this table are the condi- 
tions tested by Allen. 

The present study was undertaken to expand upon Allen’s supersonic turbulent boundary layer research by con- 
sidering a wider range of operating conditions. The results of this study, in conjunction with previous research, was 
to be used to determine what factors influence pitot probe centerline displacement in a supersonic, turbulent bound- 
ary layer. The conclusions from this paper will provide valuable insight for researchers involved in supersonic turbu- 
lent boundary layer experimentation. 
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Chapter II 

Description of Experiment 

2.1 FACILITIES 

Research was conducted in three separate NASA Lewis Research Center (LeRC) Supersonic Wind Tunnels: the 
10 x 10 ft SWT, the 15 x 15 cm SWT, and the 1 x 1 ft SWT. Details about each facility and their respective tests are 
given in the following sections. For all facilities, tunnel freestream conditions were measured using combination 
rakes in the tunnels’ plenum. 


2.1.1 Description of LeRC 10 x 10 ft SWT 

The NASA LeRC 10 ft x 10 ft Supersonic Wind Tunnel (SWT) is a propulsion tunnel primarily used to conduct 
supersonic inlet research. The 10 x 10 ft SWT is capable of producing test section Mach numbers from = 2.0 to 
M m = 3.5 with Reynolds numbers ranging from approximately 0.98xl0 6 to 11.5xl0 6 /m (0.3xl0 6 to 3.5xl0 6 /ft). 
The test section is 3.048 m high x 3.048 m wide x 12.19 m long (10x10x40 ft) and starts approximately 16.7 m (55 
ft) downstream of the throat. Figure 2.1 shows a planform view of the entire 10x10 ft SWT loop and figure 2.2 
shows an elevation view of the test section. 

2.1.2 Description of LeRC 15 x 15 cm SWT 

The 15 x 15 cm SWT tunnel provides the capability for the development and calibration of instruments at super- 
sonic conditions and for other studies related to high speed research. It has a test section 15 by 15 cm ( approxi- 
mately 0.5 x 0.5 ft). The facility is an open loop, continuous flow wind tunnel with Mach number variations 
provided by replaceable fixed-geometry nozzle blocks. Current available Mach numbers are M^ = 2.0, 2.5 and 3.0, 
and the Reynolds number can be varied from 6.0xl0 6 /m to 20xl0 6 /m. Figure 2.3 shows a sketch of the 15 x 15 cm 
SWT. 


2.1.3 Description of LeRC 1 x 1 ft SWT 

The 1 x 1 ft SWT is a continuous flow aerodynamic facility used for inlet testing, shock/boundary layer interac- 
tions, and code verification of computational fluid dynamics. The test section measures 0.305 m high by 0.310 m 
wide by 0.66 m long (1 x 1.17 x 2.17 ft). The 1 x 1 ft SWT is capable of producing Mach numbers equal to 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0, 5.0, 5.5, and 6.0. Reynolds number can be varied from 1.6xl0 6 /m to 52xl0 6 /m (0.5xl0 6 /ft to 
16xl0 6 /ft). Figure 2.4 shows a sketch of the 1 x 1 ft SWT. 


5 



2.2 TEST STATION 


2.2.1 10 x 10 ft SWT Test Station 

The 10 x 10 ft SWT test section has two sets of sidewall window inserts — one set 1.588 m (5.208 ft) from the 
beginning of the test section and the other set 4.763 m (15.208 ft) from the beginning of the test section. The 
centerlines of these windows are 254 mm (10 in.) above the tunnel centerline, and schlieren windows are normally 
mounted in these inserts. Figure 2.2 shows the window inserts. 

For this experiment, one 10 x 10 ft SWT window insert was modified to allow the insertion of a test plate which 
contained a diamond shaped hole, twelve static taps and two thermocouple plugs (fig. 2.5). The test plate was held in 
place by a retaining ring and which contained two O-ring seals to ensure proper sealing (fig. 2.6). 

The probe-tip survey station was located approximately 17.6 m (57.8 ft) downstream of the throat. The boundary 
layer at this point was turbulent and was 205 mm (8.06 in.) and 224 mm (8.83 in.) thick for Mach = 2.0 and 
Mach = 2.5, respectively. 

Outside the test section an aluminum table was mounted to the test plate frame and the test hardware was mounted 
to this table. Figure 2.7 shows a photograph of the window insert, the schlieren window, the plate insert, and the test 
table. Figure 2.8 shows a close-up of the test table and the test hardware for the 10 x 10 ft SWT test. 

2.2.2 15 x 15 cm SWT Test Station 

The 15 x 15 cm SWT test station consisted of a removable side plate that had a through hole for the probe and a 
static tap to read the local static pressure. The probe-tip survey station was located approximately 870 mm down- 
stream of the tunnel throat. The boundary layer at that point was turbulent and was equal to 1 30 mm. 

2.2 3 1 x 1 ft SWT Test Station 

The 1 x 1 ft SWT test station also consisted of a removable side plate that had a through hole for the probe and 
static taps for local static pressure. The probe-tip survey station was located approximately 195 mm downstream of 
the tunnel throat. The boundary layer at the test station varied from 25.5 to 29.3 mm, depending on the Mach num- 
ber and Reynolds number that was set. Refer back to table II for boundary layer values at each condition. 


23 INSTRUMENTATION 
23.1 10 x 10 ft SWT Instrumentation 

The 10 x 10 ft SWT boundary layer was surveyed with five round pitot probes. Four relatively large diameter 
probes — 0.32, 0.64, 0.95, and 1.27 cm (0.1875, 0.25, 0.375, and 0.5 in. respectively) — and one small diameter 
probe, 0.051 cm (0.02), were used for the survey. A photograph of the probes is shown in figures 2.9, 2.10 and 2.1 1 . 
The tips of all probes were beveled for yaw insensitivity, 2 and the average ratio of inside diameter to outside 
diameter was 0.6. The bases of the four largest probes were threaded and had O-rings for sealing. These probes 
were threaded into a separate adapter piece that was mounted in the strut (fig. 2.12). Figure 2.13 shows the 1.27 cm 
probe installed in the strut. 

The fifth and smallest diameter probe was brazed to a separate adapter and bent toward the wall to ensure that the 
probe tip touched the wall before the rest of the probe. Figure 2.14 shows the 0.051 cm diameter probe installed in 
the strut. 

The strut, shown in figure 2.15, was 660.4 mm (26 in.) long and was fabricated from stainless steel. 406.4 mm 
(16 in.) of the strut had a diamond cross sectional area and the remaining 254 mm (10 in.) had a rectangular cross 
section. A 6.35 mm (0.25 in.) diameter hole was drilled from the tip of the probe through to a slot in the rectangular 
section, and the top of the strut was milled out to allow the mounting of the adapter piece. 


2 It has been well documented (Chue, 1972) that a pitot probe with a 15° bevel is insenitive to yaw angles of ±15°. 
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The adapter piece had a 3.18 mm O.D. tube (0.125 in.) brazed to its base to produce a hard-lined route from the 
pitot tube through the strut. The adapter was fastened to the strut with a #4 screw and RTV silicon epoxy was used 
to seal around the adapter/strut interface to prevent leakage. 

The base of the pitot probe strut was bolted to a linear positioner table outside the test section. The table had a 
range of 304.8 mm (12 in.) and an accuracy of ±0.25 mm (±0.010 in.), or ±0.08 percent of full scale. The table was 
driven by a compumotor and the position of the table was read using a linear pot, and with an encoder. 

The clearance between the strut and the diamond-shaped hole in the tunnel sidewall was 0.25 mm (0.010 in.) 
total. To prevent leakage and to aide in alignment, a split-half brass sleeve was bolted to the sidewall test plate 
around the strut. The sleeve-to-sidewall plate interface had an O-ring and the sleeve-to-strut interface had O-rings. 
Vacuum grease applied to the strut helped the strut move easily through the sleeve. 

Pitot-probe pressure was read using the 10 x 10 ft SWT facility Escort D + system, a minicomputer-based, real- 
time, data acquisition system. The accuracy of this system for a 103.42 ka (15 psi) module is 0.07 percent of full 
scale. Static pressures were measured using Setra 370’ s which are more accurate at low pressures (average error of 
0.048 ka or 0.007 psia). All data were stored on a mainframe computer in the NASA-Lewis’s Research Analysis 
Center for reduction purposes. 

23.2 15 x 15 cm SWT and 1 x 1 ft SWT Instrumentation 

Both the 15 x 15 cm SWT boundary layer and the 1 x 1 ft SWT boundary layer were surveyed with five round 
pitot probes. Four “experimental” probes — 0.051, 0.091, 0.147, and 0.210 cm (0.02, 0.036, 0.058, and 0.083 in. 
respectively) — and one boundary layer survey probe, 0.036 cm (0.014). A photograph of one pitot probe is shown in 
figure 2.16. 

For both facilities, the pitot probe support tube was inserted through a grommeted hole in the sidewall then 
clamped into a Newport actuator. The grommet allowed easy movement while preventing leakage. 

Pitot probe pressures were read using the facility Escort D + system, the same type of system that exists in the 
10 x 10 ft SWT facility. All data were stored on a mainframe computer in the NASA LeRC Research Analysis Cen- 
ter for reduction purposes. 


23.3 Probe Position Zeroing Technique 

In the 10 x 10 ft SWT the pitot probe zero position was established by running the probe strut to the wall and 
using shim gages to verify that the probe tip was touching the wall. During testing, the probe was moved in 1.27 mm 
(0.05 in.) increments until the total pressure read by the probe changed (to account for any “spring” in the probe). 
The encoder value at this point was recorded by the computer and all subsequent values had this “zero” position 
value subtracted from them. 

In the 15 x 15 cm and the 1 x 1 ft SWTs the probe struts were retracted until an audible alarm sounded — this 
alarm was triggered by the stainless steel probe touching the stainless steel tunnel wall and completing an electrical 
circuit. Visual inspection (and a bit of probe manipulation) verified that the probe tip was slightly bent and thus 
touched first. With the probe tip on the wall, the computer program responsible for probe positioning was set to zero 
and the probe was actuated a couple of times to verify that the zero position was repeatable. 


2.4 Test and Procedures 

The tests were conducted in three separate NASA Lewis Research Supersonic Wind Tunnels: the 10 x 10 ft SWT, 
the 15 x 15 cm SWT, and the 1 x 1 ft SWT. Freestream conditions for each facility are given in table III. Prior to 
testing, the pitot tubes were leak checked and the zero positions of the probes were calibrated — the probes were fully 
retracted and visual inspection verified that the tip was on the wall. In the 15x15 cm SWT and 1 x 1 ft SWT facili- 
ties, the probe strut was given a small electrical charge such that when the probe tip touched the wall of the tunnel a 
circuit was completed and an audible alarm sounded. 

In the 10 x 10 ft SWT, the linear positioning table was controlled using a fluke controller— all positions were fed 
into the computer prior to testing, and the controller stepped through the program when the “Enter” key was pressed 
on the keyboard. After each position change, the static pressures beneath the probe tip were allowed to settle before 
the data point was recorded to Escort (the pitot pressure settled much faster than the static pressures). 
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In the 15 x 15 cm SWT and the 1 x 1 ft SWT the probe positions were controlled using a PACS system and a 
Newport actuator, where, again, the positions were preprogramed and the whole matrix was automated. Through 
trial and error, ample time was allowed at each y-position for the pressures to settle out before the data point was 
recorded. 
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Chapter III 

Data Reduction 

All values recorded during testing were loaded into a PC based program called PV-WAVE Personal Edition, cre- 
ated by Visual Numerics. Multiple programs were written to reduce and plot the data using the methods detailed 
below. 


3.1 EDGE MACH NUMBER 

The boundary layer edge Mach number for all the facilities was calculated from the implicit isentropic relation 
(eq. (8) of Bryer and Pankhurst, 1971): 


Pt 2 ,e 

Pt,o 


6Mg ] 

7/2 

6 

(Mg + 5 y 


7Mg — 1 


Y /2 

) 


a) 


where p t2 e is the pitot pressure measured in the freestream of the wind tunnel and p t 0 is the wind tunnel stagnation 
pressure. By calculating the Mach number by this equation, it is assumed that there is a negligible loss in stagnation 
pressure in the wind tunnel core flow between the plenum and the measurement station. This method for determin- 
ing the edge Mach number is an alternative to the more common procedure of using the freesteam pitot and local 
wall static pressures to calculate edge Mach number via the Rayleigh Pitot tube equation. The reason for using 
equation (1) was the 15 x 15 cm SWT facility did not have wall static pressure taps at the measurement station. For 
consistency, the data for the two other facilities were reduced in the same manner. 


3.2 STATIC PRESSURE 

For the purposes of data reduction, the static pressure was assumed to be constant across the boundary layer and 
was calculated from the isentropic relation: 


Ps - Pt,o 



2 

e 



( 2 ) 


For the two facilities (the 10 x 10 ft SWT and the 1 x 1 ft SWT) which had wall static pressure taps at the measure- 
ment station, the calculated and measured values were compared and found to agree within ±0.5 percent. 
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3.3 LOCAL MACH NUMBER 


Two formulas were used to calculate the local Mach number using the static pressure and the pitot pressure: 


Ps 



y-i 
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M 2 



for pressure ratios less than 1 .893 (Mach < 1 .0) and 


( 3 ) 
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for pressure ratios greater than 1 .893 (M > 1 .0). Equation (4) is the Rayleigh pitot formula from Breyer and 
Pankhurst (1971). 

The boundary layer edge Mach numbers, M e , were chosen to be the freestream Mach number. 


3.4 LOCAL VELOCITY 

Knowing the local Mach number and the edge Mach number, the local velocity ratios were then calculated from 
the following formula: 
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( 6 ) 


and, using a recovery factor of 0.89 and assuming adiabatic conditions: 

T S W „ = 0.89 (T te - T se ) + Tjg 
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The static pressure was assumed to be constant across the boundary layer, so the density ratio can be written as: 


p _ Te _l+^(0.89)M;f 

Pe T 1 + (0.89)M 1 2 


( 10 ) 


3.5 BOUNDARY LAYER TEMPERATURE DISTRIBUTION 

The temperature distribution across the boundary layer was calculated using the Crocco relation, which is valid 
for turbulent boundary layers in ideal gases. The Crocco relation is defined as: 




U 9 

1 — + 0.2M e 
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u e v u e J 


( 11 ) 


3.6 GENERALIZED VELOCITY 

The generalized velocity is the nondimensionalized velocity used for the law-of-the-wall equations. It is defined 
as such: 


U + 



where 


u* = (u e /m ,/2 jsin ^ni 1/2 u/u e j 

1 [(y-Q/^Me _ 0.2M 2 

1 1 + [(y-1)/2]m 2 1+0.2M 2 


Remember that: 
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Chapter IV 

Results 

The plots in this section are sample plots. For a complete set of plots refer to Grosser (1996). 

4.1 VERIFYING A TURBULENT BOUNDARY LAYER 

The boundary layers under investigation were verified to be fully turbulent by comparing the measured profiles 
with the theoretical law -of- the- wall for a fully turbulent boundary layer. The velocity profiles measured with the 
smallest diameter probe in each facility were transformed to law-of-the-wall coordinates using equations (12) to 
(15). The friction velocity (u T ) used to nondimensionalize the law-of-the-wall variables was deduced from Preston 
tube data obtained in each facility. The calibration procedure for determining the friction velocity is described in 
appendix B. 

An example of measured velocity profiles in law-of-the-wall coordinates is shown in figures 4.1. The theoretical 
law-of-the-wall profile using typical values for the Von Karman constant (k) and the intercept constant (C) is shown 
for comparison. The plots matches well, verifying that the boundary layers in the three NASA facilities are fully 
turbulent. 

4.2 TRUE MACH PROFILE 

To calculate the true Mach number, M Q (where here the “o” subscript denotes conditions for D = 0 probe), in each 
facility the Mach numbers of all data probes were plotted against probe diameter for each Mach number top plot of 
figure 4.2. A straight line was then fitted through the plotted points and the line was extrapolated to D = 0. The re- 
sulting M 0 points were plotted against the distance from the wall and compared to each facility’s smallest diameter 
boundary layer survey probe, (right hand plots of figs. 4.2 bottom plot in figure 4.2. For all three test sets, the com- 
parisons were good, so the calculated M o values were assumed to be the true Mach profiles. 

43 BOUNDARY LAYER THICKNESS 

In order to determine the true boundary layer thickness for each operating condition, the true Mach number pro- 
file (M 0 ) was first converted to a velocity profile (u Q ) assuming a temperature distribution based on the Crocco rela- 
tion (eq. (11)). The true boundary layer thicknesses, 8’s, were then defined as the y-positions where u/u e = 0.995. 
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4.4 CENTERLINE OFFSET 


Once the 6’s had been calculated, they were used to nondimensionalize the linear distances and all y’s were 
replaced by y/S. The diameter offset of each probe could now be shown graphically by plotting M/M e against y/5, 
as is done in figure 4.3 (fig. 4.3 is a plot of the 10 x 10 ft SWT data which exhibited the most extreme offsets). To 
calculate the actual offset values a horizontal line was drawn from each point to the ideal Mach curve figures 4.4. 
The magnitude of this horizontal line was A/5, and A was calculated simply by multiplying A/8 by 5. All of the A’s 
were plotted against y/d and an average value was determined by fitting a straight line to the majority of the points, 
figure 4.5. 

In keeping with Allen’s methods, the A for D = 0 was calculated using a different method. Allen used the slopes 
of the lines of the top graph in figure 4.2 (3M/3D) and the slopes of the true Mach number profiles in the bottom 
graph of figure 4.2 (3M/3y) and the following formulae. 

By definition, 


M - M 0 3M _ _ 

2-=> ,as D => 0 

y - y 0 d y 


but. 


y-y 0 = A 


Now, referring to Figure 4.2 
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but D 0 = 0, so: 
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plugging equations (17) and (19) into (16) leads to: 


(19) 


3D _ 9M 
A 3y 


or 


A dM/dD 
D ~ dM/dy 


( 20 ) 


at D = 0. Figure 4.6 show the results of these calculations. The points were averaged to get the A/D value for D = 0. 

A graph plotting the present study against previous studies is shown in figure 4.7. Figures 4.8 to 4.12 show the 
centerline offset error, A/D, versus D/8, M^, R TO , 3M/3y, and 8. Appendix A contains the resulting M/M e values for 
each y/8 for the three NASA facilities. The values for D = 0 are also included for comparison. 


4.5 STATISTICAL ANALYSIS 

The resulting A/D values and their corresponding Mach numbers, Reynolds numbers, and D/8 values were entered 
into RSI, a statistical analysis program by BBN Software Products, to determine which properties significantly 
influenced A/D. Table IV lists the model terms entered into RSI . 
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All properties in table IV were first normalized by RSI so that they ranged between ±1 . RSI then created more 
terms through combinations and second order terms. From the analysis RSI gave a significance for each factor, indi- 
cating which factors were not important and, thus, which factors should be removed for the next analysis. The analy- 
sis was originally run without the 15 x 15 cm SWT, R^ - 6.00x1 0 6 /m dataset (this dataset was run later) and the 
results yielded a strong freestream Reynolds number dependence and a slight freestream Mach number dependence. 
RSI generated the coefficients for a second order function, and the resulting function was: 

^ = 0.582205 + 0.057790(ML ) + 0.928536(RL ) + 0.099388(ML X R ~) + 1 78375 R L 2 (21) 

where M'^ and R'^ are the normalized freestream Mach number and normalized freestream Reynolds number and 
are defined as: 


l/R^ - 1.46681 8 x 10~ 7 
8.751022 xlO -8 


( 22 ) 


ML 


Mqq - 2.26 

0.3 


(23) 


Additional data were then taken in the 15 x 15 cm SWT tunnel at a new Reynolds number (6.00x1 0 6 /m) to see if 
its offset values conformed to the values predicted by equation (21). The new offset values did not equal the values 
predicted by equation (21), so the statistical analysis was run again, this time with the new 15 x 15 cm SWT data 
included. Instead of freestream Reynolds number, Reynolds number based on probe diameter, R D , was used (actu- 
ally, R^ was first used, then switched to R D ). RSI indicated that R D , D/5, and 5 were not significant factors in the 
pitot probe centerline displacement. It indicated that with the data given, the significant factors were Mach gradient, 
t)M/3y, and freestream Mach number, M^. The calculated coefficients yielded the following equation: 


— = 1.453183 + 0.1 87443(ML ) + 0.9981 — ) - 0.748 1 58| — 
d l J V oy 


where R’^ and M’^ are defined in equations (22) and (23) and 5M'/dy is defined as: 


(3M70y) = 


(l/(dM/3y))- 1.3094945 
0.01195055 


same as before. Figure 4. 1 3 shows a three-dimensional plot of these results. 


(24) 


(25) 
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Chapter V 

Conclusion 


The purpose of this research was to create a supersonic, turbulent boundary layer dataset and to attempt to define 
the factors influencing pitot probe centerline displacement error in such a flowfield. An additional purpose was to 
formulate the relationship between the offset error and its factors. This research has indeed created a dataset for this 
flowfield, but it has also indicated the need for more data. The data taken proved to be a small part in a very large 
and complex picture. Nonetheless, the statistical analysis program RSI was used to fit the existing data and gener- 
ated the following formula: 


A. 

D 


1 .453 1 83 + 0. 1 87443(NC ) + 0.998 


( m ') 

-0.748158 

(dM') 

J 


l ay J 


where the primed values are defined in equations (22) and (25). 

Figures 4.08 to 4.12 plot all the factors separately and figure 4.13 gives a three-dimensional graph of this formula. 
Referring to figure 4.1 1 it looks as if the error relation increases asymptotically as dM/dy goes to zero, but logically, 
that is not true. In reality, a point must exist at 0.0. 0.0 indicating that in the freestream the probes read the correct 
pressures (which they did in all experiments). This generates some interesting questions about what is happening to 
the correction curve between dM/dy = 0.000/mm and 3M/3y = 0.030/mm. Unfortunately, the statistical analysis 
program RSI could not generate an accurate function with the 0.0, 0.0 point without more data. Figure 5.1 illustrates 
the possible curves RSI might generate without additional data to fill the 0.000/mm < 3M/0y < 0.030/mm gap. An- 
swers can only be found with more experimentation. 

Another noteworthy point is that as Mach number decreases, and, specifically, when M<1.0, A/D becomes less 
than one. This is very important because it indicates that the shock in the supersonic boundary layer is skewing the 
flow and causing the probe to read a pressure outside the diameter of the probe (something that did not seem physi- 
cally possible). Again, more data is needed to define exact limits. 

The results from this research can be summarized as thus: 

1 . In a supersonic turbulent boundary layer the pitot probe centerline offset error, A/D, is not greatly influenced 
by the probe diameter, D/8, the boundary layer thickness, 5, or the Reynold numbers R^ or R D 

2. The pitot probe centerline displacement error is influenced by the Mach gradient, 3M/3y, and the freestream 
Mach number, M^. 

3. A/D = 1 .453183 + 0.187443 (M 1 ^) + 0.998(3MY3y) - 0.748 158(3M73y) 2 3 4 appears to describe the relationship 
between offset error and test conditions for 0.030/mm <dM/dy <0.1 00/mm. Primed variables are defined in equa- 
tions (22) and (25). This equation must go to zero as 3M/0y goes to zero. 

4. More data is needed to complete the supersonic, turbulent boundary layer matrix and to allow for more accu- 
rate formulation. 
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5. Future work should include a repeat of the 10x10 ft SWT = 2.00 and M w = 2.56 data and should also at- 
tempt to fill the gaps between 0.000/mm <3M/3y <0.030/mm. This may be achieved by running the 10 x 10 ft SWT 
at higher Mach numbers. 

6. It appears that for most supersonic, turbulent boundary layer research the correction factor lies near zero. Thus, 
for research in tunnels where 5M/3y is greater than 0.030/mm, a correction factor is not needed. 
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APPENDIX A 


SUMMARY OF REDUCED DATA 


Summary of Reduced 10 x 10 ft SWT Data; M m - 2.00, = 4.27x 1 Ofym 



M/Me 


y/delta 

D=12.7 cm 

D=9.525 cm 

D=635 cm 

D=3.175 cm 

D=0 cm 

0.06 

■gglgra 

0.667 


0.633 

0.623 

0.11 

■ 

0.717 

0.706 

0.690 

0.678 

0.17 

0.762 

0.751 

0.742 

0.726 

0.716 

0.22 

0.789 


0.770 

0.755 

0.745 

0.28 

0.815 

1' 

0.795 

0.781 

0.771 

0.34 

0.839 


0.821 

0.804 

0.795 

0.39 

0.862 


0.841 

0.827 

0.816 

0.45 

0.881 


0.862 

0.848 

0.838 

0.51 

0.900 


0.883 

0.868 

0.860 

0.56 

0.918 


0.901 

0.887 

0.879 

0.62 

0.933 


0.917 

0.906 

0.897 

0.67 

0.949 

0.940 

0.934 

0.924 

0.917 

0.73 

0.961 

0.953 

0.949 

0.940 

0.934 

0.79 

0.972 

0.966 

0.961 

0.955 

0.950 

0.84 

0.981 

0.977 

0.974 

0.967 

0.963 

0.90 

0.989 

0.986 

0.984 

0.980 

0.977 

0.96 

0.993 

0.991 

0.990 

0.988 

0.986 

1.00 

0.996 

0.995 

0.995 

0.992 

0.991 

1.01 


0.996 

0.996 

IBXEEHI 

■i MM 

mm 

0.999 

0.998 

1.000 


— 

Ira 

1.000 

1.000 

1.001 


0.999 

Ira 

1.001 

1.001 

1.002 

1.001 

1.002 

Ira 

1.002 

1.002 

1.002 

1.002 

1.003 

m 

1.002 

1.003 

1.003 

1.003 

1.004 

1.35 

1.003 

1.003 

1.003 

1.003 

1.003 

1.41 

1.003 

1.003 

1.003 

1.003 

1.003 
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Summary of Reduced 10 x 10 ft SWT Data; M B = 2.56, = 4.27x1 0 6 /m 



M/Me j 


y/delta 

D=12.7 cm 


D=6.35 cm 

D=3.175 cm 

D=0 cm | 

0.04 

0.609 


0.574 

mssm 

MSSSM 

0.08 

0.666 

0.642 

0.629 



0.12 

0.693 

0.674 

0.664 

0.645 


0.17 

0.719 

0.701 

0.690 

0.674 

0.660 

0.21 

0.744 

0.724 

0.715 

0.699 

0.685 

0.25 

0.766 

0.746 

0.737 

0.721 

0.706 

0.29 

0.787 

0.767 

0.756 

0.743 

0.727 

0.33 

0.807 

0.787 

0.776 

0.761 

0.746 

0.37 

0.828 

0.807 

0.797 

0.781 

0.765 

0.41 

0.847 

0.825 

0.816 

0.802 

0.786 

0.46 

0.865 

0.843 

0.833 

0.821 

0.805 

0.50 

0.883 

0.861 

0.851 

0.839 

0.823 

0.54 

0.898 

0.878 

0.868 

0.857 

0.841 

0.58 

0.914 

0.894 

0.885 

0.873 


0.62 

0.929 

0.909 

0.901 

0.890 


0.66 

0.945 

0.924 

0.917 

0.906 

0.892 

0.70 

0.958 

0.939 

0.932 

0.922 

0.909 

0.74 

0.967 

0.952 

0.944 

0.937 

0.925 

0.79 

0.977 

0.963 

0.957 

0.950 

0.940 

0.83 

0.984 

0.974 

0.967 

0.962 

0.954 

0.87 

0.989 

0.982 

0.975 

0.971 

0.964 

0.91 

0.993 

0.987 

0.984 

0.981 

0.976 

0.95 

0.996 

0.992 

0.990 

0.987 

0.984 

0.99 

0.998 

0.995 

0.993 

0.991 

0.988 

1.00 

0.998 

0.996 

0.993 

0.991 

0.989 

1.03 

0.999 

0.998 

0.996 

0.995 

0.993 
























Summary of Reduced 1 x 1 ft SWT Data; - 1.96, R„ = 16.9xl0 6 /m 
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Summary of Reduced 1 X 1 ft SWT Data; = 1.96, R„ = I0.3xl0 6 /m 



M/Me | 


y/delta 

D =0.2 10 cm 

D=0.147 cm 

D =0.090 cm 

D=0.051 cm 

D=0 cm 




0.609 


0.607 




0.649 


0.646 

0.17 


1 

0.681 

■ 

0.679 


0.711 

0.710 

0.711 

1 

0.709 


0.738 

0.735 

0.733 

0.734 

0.731 

0.30 

0.758 

0.757 

0.756 

0.757 

0.756 

0.35 

0.782 

0.780 

0.778 

0.778 

0.777 

0.39 

0.803 

0.802 

0.800 

0.800 

0.799 

0.44 

0.825 

0.823 

0.820 

0.821 

0.819 

0.48 

0.844 

0.844 

0.842 

0.842 

0.841 

0.53 

0.864 

0.864 

0.863 

0.863 

0.862 

0.57 

0.885 

0.884 

0.884 

0.883 

0.883 

0.62 

0.903 

0.901 

0.903 

0.902 

0.902 

0.66 

0.920 

0.921 

0.921 

0.919 

0.919 

0.71 

0.937 

0.938 

0.938 

0.936 

0.937 

0.75 

0.953 

0.952 

0.953 

0.952 

0.952 

0.80 

0.966 

0.965 

0.966 

0.965 

0.965 

0.84 

0.976 

0.975 

0.976 

0.975 

0.975 

0.89 

0.984 

0.984 

0.983 

0.982 

0.983 

0.93 

0.989 

0.989 

0.989 

0.988 

0.988 

0.98 

0.992 

0.992 

0.992 

0.991 

0.991 

1.00 _ 

0.993 

0.993 

0.993 

0.992 

0.992 





0.993 

HU 

1 



■ 

0.994 

■S 

1.11 

1 


■ H 

0.995 

mm 

1.16 

| 0.997 

0.996 

0.996 

0.996 

■H 


Summary of Reduced 1 x 1 ft SWT Data; = 2.46, R„ = 10.3xl0tym 



M/Me | 


y/delta 

D=0J10 cm 

D=0.147 cm 


D=0.051 cm 

D=0 cm 



0.567 



0.565 



0.607 




0.18 


0.641 

0.640 

0.645 

0.641 


0.672 


0.666 

0.673 

0.669 


0.696 

0.693 

0.690 

0.697 

0.693 

0.32 

0.721 

0.719 

0.714 

0.721 

0.717 

0.37 

0.746 

0.743 

0.738 

0.745 

0.741 

0.42 

0.768 

0.767 

0.763 

0.769 

0.766 

0.47 

0.793 

0.790 

0.787 

0.792 

0.789 

0.52 

0.815 

0.814 

0.810 

0.815 

0.812 

0.56 

0.836 

0.835 

0.834 

0.837 

0.836 

0.61 

0.859 


0.855 

0.859 

0.857 

0.66 

0.880 


0.877 

0.880 

0.878 

0.71 

0.900 

0.899 

0.899 

0.900 

0.899 

0.76 

0.921 

0.918 

0.918 

0.920 

0.918 

0.80 

0.939 

0.937 

0.937 

0.938 

0.937 

0.85 

0.954 

0.953 

0.953 

0.955 

0.955 

0.90 

0.969 

0.968 

0.969 

0.970 

0.970 

0.95 

0.981 

0.979 

0.980 

0.981 

0.981 

1.00 

0.989 

0.988 

0.989 

0.989 

0.989 

1.00 

0.989 

0.988 

0.989 

0.990 

0.990 

1.04 

0.994 



0.995 


1.09 

0.996 



0.997 


1.14 

0.998 

■ 

0.998 

0.999 

0.999 

1.19 

0.998 

1 

0.999 

0.999 

1.000 

1.24 

0.998 

: . ■ - 

0.999 

0.999 

1.000 











































Summary of Reduced 15 x 15 cm SWT Data; M m = 2.47, R w = 10.3xl0 6 /m 



M/Me 


y/delta 

0=0.210 cm 

D=0.147 cm 

D=0.090 cm 

1 

0 

1 

D=0 cm 

0.12 

0.592 

0.589 

0.591 

0.593 

0.592 

0.17 

0.638 

0.624 

0.625 

0.627 

0.620 

0.21 

0.669 

0.655 

0.654 

0.655 

0.648 

0.26 

0.689 

0.679 

0.679 

0.680 

0.675 

0.30 

0.714 

0.699 

0.703 

0.706 

0.700 

0.35 

0.732 


0.727 

0.729 

0.726 

0.39 

0.753 

0.749 

0.752 

0.750 

0.749 

0.44 

0.779 

0.770 

0.776 

0.772 

0.770 

0.48 

0.802 

0.790 

0.796 

0.791 

0.788 

0.53 

0.823 

0.810 

0.816 

0.812 

0.809 

0.58 

0.842 

0.836 

0.835 

0.835 

0.832 

0.62 

0.861 

0.859 

0.859 

0.855 

0.854 

0.67 

0.882 

0.869 

0.880 

0.876 

0.875 

0.71 

0.904 

0.897 

0.896 

0.898 

0.894 

0.76 

0.921 

0.920 

0.918 

0.914 

0.913 

0.80 

0.938 

0.934 

0.938 

0.933 

0.933 

0.85 

0.953 

0.944 

0.952 

0.949 

0.949 

0.89 

0.963 

0.960 

0.970 

0.962 

0.966 

0.94 

0.973 

0.975 

0.980 

0.974 

0.977 

0.98 

0.986 

0.983 

0.990 

0.984 

0.987 

1.00 

0 991 

0.986 

0.994 

0.988 

0.990 




0.999 



1.07 

0.999 


1.002 


1 

1.12 

1.002 

0.999 

1.002 


1.002 

1.16 

1.002 


1.002 

1.001 

1.001 

1.21 

1.002 

1.001 

1.001 

1.001 

1.001 


Summary of Reduced 15 x 15 cm SWT Data; = 2.47, R^ = 6.00x1 0 6 /m 



M/Me 


y/delta 

D=0.210 cm 

D=0.147 cm 

D=0.090 cm 

D=0.051 cm 

D=0 cm 

0.12 

0 567 


0.569 

0.578 

0.578 

0.16 

0.609 


0.604 

0.610 

0.607 

0.20 

0 639 


0.629 

0.635 

0.630 

0.24 

0.669 


0.654 

0.661 

0.653 

0.28 

0.688 


0.677 

0.685 

0.679 

0.32 

0,711 


0.701 

0.706 

0.701 

0.36 

0.728 


0.723 

0.727 

0.724 

0.41 

0.750 


0.744 

0.751 

0.747 

0.45 

0.774 


0.764 

0.771 

0.765 

0.49 

0.795 


0.785 

0.792 

0.787 

0.53 

0.814 

■ 

0.806 

0.810 

0.806 

0.57 

0.833 


0.827 

0.831 

0.828 

0.61 

0.852 

0.845 

0.847 

0.853 j 

0.850 

0.65 

0.871 

0.870 

0.867 

0.871 | 

0.868 

0.69 

0.889 

0.887 

0.863 

0.889 | 

0.886 

0.73 

0.910 

0.905 

0.901 

0.907 

0.902 

0.78 

0.929 

0.921 

0.920 

0.922 

0.918 

0.82 

0.942 

0.938 

0.937 

0.940 

0.937 

0.86 

0.958 

0.952 

0.950 

0.955 

0.950 

0.90 

0.969 

0.965 

0.967 

0.968 

0.967 

0.94 

0.981 

0.977 

0.980 

0.978 

0.978 

0.98 

0.991 

0.987 

0.988 

0.987 

0.985 

1.00 

0.994 

0.991 

0.991 

0.991 

0.990 

1.02 

0.997 

0.995 




1.06 

1.003 

1.001 



■ 

1.10 

1.006 

1.005 

1.005 

1.005 

■ 










































APPENDIX B 


CALCULATION OF C f 


CALCULATION OF C f 

C f is defined as (2x w )/(pu) and was calculated by using the Bradshaw-Unsworth (1974), equation. This equation 
was chosen over Sigalla (1965) and Fenter-Stalmach (1957) because it is based on the law-of-the-wall equations and 
not on the power- law equations. Bradshaw and Uns worth used a Preston tube — a circular pitot tube resting on the 
surface — to measure the total pressure at the wall. This pressure is then plugged into an equation and the equation is 
solved for C f . The Bradshaw-Unsworth equation is as follows: 


P t =96 + 6()log l0 ||i| + 2.37 


log.ol^ 


+ 10 M: 


:( r o.3° 



where 


Pt = 


P t - p. 

wall 

^wall 


R x = 


D 

Fw 




m t 
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TABLE I.— SUMMARY OF PREVIOUS RESEARCH 


Involigator 

CVS 

Midi 

Number 

Reynokb 
Number 
per mtter 

Flow 

Candiliom 

General Reautn 

Youn**. Man 

0.02 to 0.2? 

inooiqpreHibk 

1.3 k 10* 

Wake 

ad-o.u 

F. Daviea 

0.18 to 11* 

2.3 

10.0x10* 

to 

12.0x10* 

Laminar B L 
(flat plate ft oooa) 

No 

dupUomMt 

Liveaay 

0 01 to 0.11 

InranprawMi 

0.5 x 10* 

Torbulctt B-L 

AD-0.09 

P. Davitt 

0 04 to 0.69 

Inoomprwalbte 

1.5 X 10* 

Turbulent B L 
Wake 

N..M— 

MacMJlu 

0.02 to 0.12 

inoompnanbk 

0-5 to 2.3 x 10* 

Turbulent B L 
TurbulerX Pipe 

ATM). 15 for 
y/d>2 

Alien 

0.018 to 0.691 

2.0 

8 x 10* 

Turbulent B L 

ATM). 31 for 
D-0 


TABLE II. — SUMMARY OF RESEARCH VARIABLES 


Facility 

Nom. 

Mach 

Frocatraua 

Reynold* 

Number 

Approx. 

Boundary 

Layer 

Thlckneai 

D/8 

lOftx 10 ft SWT 

2.00 

4.27 x l0‘/m 

224 mm 

.015 

.030 

.045 

.060 


256 


205 mm 





15 cm x 15 cm 

2.47 

10 3 x 10* An 

12.7 mm 

039 

.072 

.121 

.167 

SWT 


6.00 x lO'/m 

13.0 mm 





1 ft x 1 ft SWT 

1 96 

10 3 x 10* /m 

29.3 mm 

.016 

.030 

.053 

.070 


2.46 


27.5 mm 






1 96 

16.9xl0‘/m 

26.4 mm 






2 46 


25.5 mm 





4 ft x 4 ft** SWT 

20 

8.0 x 10* /m 

69.9 mm 

018 

.034 

.145 

.256 

(Allen) 









** Alien also ran D/6 - .364, .473, .582, .691 . 


TABLE III.— TUNNEL FREESTREAM CONDITIONS 


Facility 

M 

Re. (/m) 

P,(ka) 

T,(°C) 

10 ft x 10 ft - SWT 

2.00 

4.26 x 10* 

33.8 

27 


2.56 


44.8 


15 cm x 15 cm 

2.47 

10.3 x 10* 


mm 

SWT 


6.00 x 10* 


warn 

1 ft X 1 ft SWT 

1 96 

10.3 x 10‘ 




2.46 

16 9 x 10* 



4 ft x 4 ft SWT 

2.00 

8.0 x 10 4 

202.6 

43 

(Allen) 






* For the 0 05 1 cm diameter probe at M=2 .56, P, = 56.2 ka and T, = 37 °C 
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TABLE IV.— LIST OF RSI MODEL TERMS 


Item 

Model Term 

1 

D/5 

2 

M, 

3 

R- 

4 

dM/dy 

5 

5 


TABLE V.— LIST OF A/D, AND D/5 VALUES 



A/D 

D/5 

Mach 

Number 

Reynold! 
Number 
per meter 


5 

ran 


0.962 

0.000 

2.00 

4.27E+06 

0.004 

224.0 


1.321 

0.014 






1.693 

0.029 





10 ft x 10 ft Data 

1 484 

0.042 






1.594 

0.057 






1.661 

0.000 

2.56 

4.27E+06 

0.006 

205.0 


1.984 

0.016 






1.721 

0.031 






1.864 

0.046 






2 083 

0.062 






0.088 

0.000 

1.96 

10.3E+06 

0.029 

29.3 


0.033 

0.017 






0.090 

0.031 






0.064 

0.050 





1 ft x 1 ft Da La 

0.072 

0.072 






0.029 

0.000 

2.46 

10.3E+06 

0.041 

27.5 


0.220 

0.019 






-0.122 

0.033 






-0.012 

0.053 






0.049 

0.076 


WSm MM 




-0.123 

ci'ooo 

1.96 

Btrli lti ill 

0.028 

26.4 


0.108 

0.019 






-0 369 

0.034 


■ 




-0.228 

0.057 


■ 



1 ft % 1 ft Data 

-0.140 

0.080 


mmmmm 




-0.096 

0.000 

2.46 

16.9E+06 

0.040 

25.5 


-0.205 

0.020 






-0.1 14 

0.036 






-0.071 

0.058 






-0.114 

0.083 






0.062 

0.000 

2.47 

10.3E+06 


12.7 


0.101 

0.166 





15cm x 15 cm Data 

-0.006 

0.116 






0.130 

0.072 






0.083 

0.040 






0.043 

0.000 


6.00E+06 

0.091 

13.0 


0.085 

0.166 





15 cm x 15 an Data 

0.010 

0.116 






-0.021 

0.072 






0.203 

0.040 





am 

0.480 

0.000 

■FTTiS 

8.00E+06 

0.020 

699 


0.500 

0.010 






0.520 

0.034 
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u/U «o/D 


0.140 

0.120 

0.100 

0.080 

0.060 

0.040 

0.020 


0.000 


0.000 



A 


R p = Reynolds number based on pipe diameter 


0.050 0.100 


Figure 1.2. — MacMillan’s apparent displacement 8 0 . 
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A/D 
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Static taps 


■- Start of test 
' section 



Slot for sting 


Personnel 
access 
door - 



Figure 2.2. — Elevation view of the 10 ft x 10 ft SWT Test section. 
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Figure 2.3. — Sketch of NASA LeRC 1 5 cm x 1 5 cm SWT. 


1FT X 1FT VARIABLE MACH WIND TUNNEL (BUILDING 37 - ROOM 1 NW) 
MACH NUMBER 1.3 THRU 6.0 — MAXIMUM TEMPERATURE 650 *F 

r mm SUPPLY PRESSURE • 450 PS! 

va mr 

TEST SECTOR- ... MZZl E AIR INLET • 40 PH * 125 PSI 



C-89-07886 

Figure 2.4. — Sketch of NASA LeRC 1 ft x 1 ft SWT. 
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Figure 2.7. — Far view of test hardware outside oflOftxIOft SWT test section. 



Figure 2.8. — Close-up view of test hardware on support table. 
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pjfjl 


C -95-00275 


Figure 2.10. — Photograph of pitot probe leading edges. 
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C-95-00280 

Figure 2.1 1 . — Photograph of diameter = 0.051 cm pitot probe. 



Figure 2.1 2. — Photograph of adapter piece installed in research strut. 
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Figure 4.1 . — Example of law of the wall profiles. 



D, mm 
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M/M 0 M/M ( 



True profile 

O d/8 = 0.066 

(D = 1 .270 cm) 


True profile 
D/8 = 0.050 
(D = 0.952 cm) 


0.5 1.0 1.5 2.0 

y/8 




I 0.8 ^ 


True profile 


True profile 

D/8 = 0.033 

0.6 ® 

O D/8 = 0.017 

(D = 0.635 cm) 


(D = 0.317 cm) 


0.4 


1.0 1.5 2.0 

0.0 

0.5 1.0 1.5 2 


Figure 4.4. — Example of calculation of pitot probe centerline offset errors. 
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o 


A/D ave = 0.043 


Mach = 2.47 


(for D =0) 


0.2 0.4 0.6 

y/8 

Figure 4.6.— Example of calculation of A at diameter = 0. 


0.8 


O M = 2.00, R e = 4.27 x 1 0^ /m, 8 = 22.4 cm 
♦ M = 2.56, R e = 4.27 x 1 06 /m, 8 = 20.5 cm 




O M = 2.00 (Allen), R e = 8x10 6 /m, 8 = 6.99 cm 
■ M = 2.47, Re = 1 .03 x 1 0 7 /m, 8 = 1 .27 cm 

▲ M = 2.47, Re = 6.00x1 06 /m, 8 = 1.30 cm 

v M = 1 .96, Re = 1 -03 x 1 0 7 /m, 8 = 2.93 cm 

▼ M = 2.46, R e = 1.03x1 0 7 /m, 8 = 2.75 cm 

o M = 1.95, R e = 1.69 x 10 7 /m, 8 = 2.64cm 

• M = 2.46, R e = 1 -69 x 1 0 7 /m, 8 = 2.55 cm 


0.02 0.04 0.06 0.08 


D/8 

Figure 4.7. — Effect of probe size on displacement ratio in boundary layers. 



A/D 




o 

o 

o 



-0.500 L 

0 2 4 6 8 10 12 14 16 18: 

Freest ream Reynolds number, R^, /m 


Figure 4.10. — Effect of freestream Reynolds number, R^, on centerline offset error. 


O GOO O 






^ 1.000 


o 




Figure 4.13. — 3-Dimentional graph of A/D versus FL and 



e 


2.500 
2.000 

1.500 
1.000 
0.500 
0.000 

-0.500 



0.000 0.010 0.020 0.030 0.040 0.050 0.060 0.070 0.080 0.090 0.100 


Mach gradient, 2M/2y, /mm 


Figure 5.1. — Possible correction curves. 


46 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 


1. AGENCY USE ONLY (Leave blank) 


4. TrTLE AND SUBTITLE 


2. REPORT DATE 


January 1997 


REPORT TYPE AND DATES COVERED 

Technical Memorandum 


Factors Influencing Pitot Probe Centerline Displacement in a Turbulent 
Supersonic Boundary Layer 


6. AUTHOR(S) 

Wendy I. Grosser 


FUNDING NUMBERS 


WU- 505-62-52 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


B. PERFORMING ORGANIZATION 
REPORT NUMBER 


E- 10484 


9. SPONSORING/MON ITO RING AGENCY NAME(S) AND ADDRESSEES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 107341 


11. SUPPLEMENTARY NOTES 


This report was submitted as a dissertation in partial fulfillment of the requirements for the degree Bachelor of Science in 
Aerospace Engineering to the University of Minnesota in Minneapolis, Minneapolis, Minnesota 55455-01 10. Responsible 
person, Wendy I. Grosser, organization code 2830, (216) 433-8680. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 35 

This publication is available from the NASA Center for AeroSpace Information, (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

When a total pressure probe is used for measuring flows with transverse total pressure gradients, a displacement of the 
effective center of the probe is observed (designated A). While this phenomenon is well documented in incompressible 
flow and supersonic laminar flow, there is insufficient information concerning supersonic turbulent flow. In this study, 
three NASA Lewis Research Center Supersonic Wind Tunnels (SWT’s) were used to investigate pitot probe centerline 
displacement in supersonic turbulent boundary layers. The relationship between test conditions and pitot probe centerline 
displacement error was to be determined. For this investigation, ten circular probes with diameter-to- boundary layer ratios 
(D/8) ranging from 0.015 to 0.256 were tested in the 10 ft x 10 ft SWT, the 15 cm x 15 cm SWT, and the 1 ft x 1 ft SWT. 
Reynolds numbers of 4.27 x 10 6 /m, 6.00 x 10 6 /m, 10.33 x 10 6 /m, and 16.9 x 10 6 /m were tested at nominal Mach 
numbers of 2.0 and 2.5. Boundary layer thicknesses for the three tunnels were approximately 200 mm, 1 3 mm, and 30 
mm, respectively. Initial results indicate that boundary layer thickness, 8, and probe diameter, D/ 8 play a minimal role in 
pitot probe centerline offset error, A/D. It appears that the Mach gradient, OM/ 8y, is an important factor, though the 
exact relationship has not yet been determined. More data is needed to fill the map before a conclusion can be drawn with 
any certainty. This research provides valuable supersonic, turbulent boundary layer data from three supersonic wind 
tunnels with three very different boundary layers. It will prove a valuable stepping stone for future research into the factors 
influencing pitot probe centerline offset error. 


14. SUBJECT TERMS 


Turbulent boundary layer; Pitot probe displacement; Supersonic 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE 

Unclassified Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

48 


16. PRICE CODE 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 




























